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In this paper, the holographic dark energy model with new infrared (IR) cut-off for both the 
flat case and the non-flat case are confronted with the combined constraints of current cos- 
mological observations: type la Supernovae, Baryon Acoustic Oscillations, current Cosmic Mi- 
crowave Background, and the observational hubble data. By utilizing the Markov Chain Monte 
Carlo (MCMC) method, we obtain the best fit values of the parameters with la, 2a errors in 
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n m0 = 0.2762+^°32olo.oli_, H = 70.2Q±f™ttf . In the best fit holographic dark energy models, 
| the equation of state of dark energy and the deceleration parameter at present are characterized 

! by w de0 = -1.1414 ± 0.0608, go = -0.7476 ± 0.0466 (flat case) and w de0 = -1.0653 ± 0.0661, qo = 

£SJ ' —0.6231 ± 0.0569 (non-flat case). Compared to the ACDM model, it is found the current combined 

datasets do not favor the holographic dark energy model over the ACDM model. 
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Since 1998, the type la supernova (SNe la) observations [l|, have shown that our universe has entered into a 
phase of accelerating expansion. During these years from that time, many additional observational results, including 
, current Cosmic Microwave Background (CMB) anisotropy measurement from Wilkinson Microwave AnisotropyProbe 
(WMAP)d, i], and the data of the Large Scale Structure (LSS) from Sloan Digital Sky Survey (SDSS) 0,11], also 
strongly support this suggestion. These observational results have greatly inspirited theorists to understand the 
mechanism of the accelerating expansion of the universe, which is usually attributed to an exotic energy component 
with negative pressure, dubbed dark energy (DE). The simplest but most natural candidate of DE is the cosmological 
constant A, with the constant equation of state (EOS) w — —1. As we know, the cosmic concordance model confronts 
with two difficulties: the fine-tuning problem and the cosmic coincidence problem. Both of these problems are related 
to the DE density. In order to solve or alleviate cosmological constant puzzles, many dynamical DE models are 
proposed, where the DE density and its EOS are time-varying. However, the predictions of the cosmological constant 
model still fit to the current observations @-[_i- Therefore the dynamical DE models being proposed should not be 
far away from the cosmological constant model, such as quintessence [T3___| > phantom [Hp j_quintom [I?} i K-essence 
[TU, tachyon [!_}, ghost condensate [__}, holographic DE [_l|, __| and agegraphic DE [__{ __| etc. Although many 
DE models have been presented, the nature of DE is still a conundrum. This puzzle can not be understood before 
a complete theory of quantum gravity is established. But the two additional aspects from the current cosmological 
observations and some basic quantum gravitational principles may shed light on probing the nature of DE. 

On the one hand, provided that we know little on the theoretical nature of DE at present, the combined cosmic 
observations can play an important role in understanding the nature of DE. The cosmological parameters space in the 
DE model can be determined by the constraints of the data combinations. Recently, the 397 SN la data was compiled 
in Ref. [__} by adding CfA3 sample from the CfA SN Group to the Union set by Ref. [__}, which include 250 SN la 
at high redshift but only 57 at low redshift, to form the Constitution set. Aside from the SN la data, the combined 
analysis is required in order to break the degeneracy between the cosmological parameters, which includes cosmic 
observations from baryon acoustic oscillations (BAO), CMB and the observational Hubble data (OHD). The BAO are 
detected in the clustering of the combined 2dFGRS and SDSS main galaxy samples or the SDSS luminous red galaxies 
and measure the distance-redshift relation. From these samples, the values of [r s (zd)/ Dv{0.2), r s (zd)/ Dy (0.35)] and 
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their inverse covariance matrix in the measurement of BAO can be obtained [27j. For the measurement of CMB, 
we utilize the shift parameter R at the photon decoupling epoch z*, the acoustic scale Ia( z *), and together with the 
physical baryon density parameter multiplied by 100, thus it is lOOfl^h 2 [H, [2§|. Here, it is worth noting that the 
WMAP distance information R(z*) and Ia(z*) can not be measured by WMAP directly, but are derived from making 
a global fitting constraint with MCMC method by using the full WMAP data on the assumption that a certain 
cosmological model has been given in advance [30l |. Although in theory the inverse covariance matrix on R(z*) and 
Ia{z*) is model dependent, it is feasible to use the derived results about R(z*) and Ia(z*) to constrain the parameters 
in another DE model since R(z*) and Ia(z*) do not depend strongly on the DE model which is not far away from the 
cosmological constant model [3(|. What is more, the paper [3l[ has been demonstrated that [R(z*), Ia(z*), WOflbh 2 ] 
effectively provide a good summary of CMB data when the DE model parameters are constrained. In addition, we 
employ the OHD at twelve different redshifts determined by using the differential ages of pas sively evolving galaxies in 
Ref. [33], where the value of the Hubble constant is replaced by Ho = 74.2 ± 3.6 in Ref. [33|], and add the three more 
observational data H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H±z = 0.43) = 86.45 ± 3.27 in 0. 
Since the constraint results of a given model are dependent on the combined data (35l . l36l ] , in this paper we use a fully 
combined observations from the 397 SN la standard candle data, the value of [r s (zd)/ Dy (0.2), r s (zd) / Dy (0.35)) and 
their inverse covariance matrix in the measurement of BAO, the values of [R(z*), Ia(z*), lOOVl^h 2 ] and their inverse 
covariance matrix in the measurement of CMB, and the fifteen OHD. 

On the other hand, the models which are constructed in light of some fundamental principle are more charming, 
since this kind of DE model may exhibit some underlying features of DE, for instance the holographic DE model |2ll.[22| 
and the agegraphic DE model [H, [24| ■ The holographic DE model is built on the basis of holographic principle and 
some features of quantum gravity theory. The agegraphic DE model is derived from taking the combination between 
the uncertainty relation in quantum mechanics and general relativity into account. In this paper, we focus on the 
holographic DE model, which is considered as a dynamic vacuum energy. According to the holographic principle, the 
number of degrees of freedom in a bounded system should be finite and is related to the area of its boundary. By 
applying the principle to cosmology, one can obtain the upper bound of the entropy contained in the universe. For 
a system with size L and UV cut-off A without decaying into a black hole, it is required that the total energy in a 
region of size L should not exceed the mass of a black hole of the same size, thus L 3 p^ < LM%. The largest L allowed 
is the one saturating this inequality, thus we obtain the holographic DE density 

3c 2 A/ 2 , 

PA = -^, (1) 



where c is a numerical constant and M p i is the reduced Planck Mass M p i = l/y8nG. It just means a duality between 
UV cut-off and IR cut-off. The UV cut-off is related to the vacuum energy, and IR cut-off is related to the large scale 
of the universe, for example Hubble horizon, particle horizon, event horizon, Ricci scalar or the generalized functions 
of dimensionless variables as discussed by [2l|, H3, HM|- Next, we give a brief review on the main results when 
Hubble horizon, particle horizon, event horizon or Ricci scalar are taken as the IR cut-off, respectively. 

• L~ 2 = H 2 . As pointed in [22|], it is found that the holographic DE density is in proportion to H 2 , the same as 
dark matter density, i.e. pde/c 2 = Pmj(\ — c 2 ) oc H 2 . It appears that it is natural to solve the coincidence problem. 
However, Hsu [2l[ pointed out that the dark energy EOS Wde — was obtained in this instance. It is obvious that 
this result is not consistent with the current observations. This bad situation can be changed by considering the 
holographic DE with Hubble horizon as the time variable cosmological constant. More detailed analysis is presented 
in Ref. [39j |. 

• L~ 2 = Rph(a) = a f Q = a J ° j^m- As shown in paper (22[, Li pointed out that this yields the dark energy 
EOS is not less than —1/3. Thus the current accelerated expansion of our universe can not be well explained. However, 
this result in [22| is obtained on the assumption that DE dominates. The holographic DE model with particle horizon 
has been discussed in detail by [4p| . 

• L~ 2 = R e i l (a) = a -^j^ = a f a °° ■ The holographic DE model with event horizon can reveal the dynamic 
nature of the vacuum energy and provide a desired EOS of the holographic DE with the model parameter c. Further- 
more, the holographic DE behaves like quintessence, cosmological constant and phantom respectively for the different 
values of the model parameter: c > 1, c = 1 and c < 1 [22] . Therefore, the value of model parameter c plays a crucial 
role in determining the property of holographic DE in this case. However, this model is confronted with the causality 
problem: why should the present density of DE be determined by the future event horizon of the universe. 

• L- 2 = R ! = -6(H + 2H 2 + \). In [33], it has shown that this model can avoid the causality problem and 
naturally solve the coincidence problem of dark energy after Ricci scalar is taken as the IR cut-off and the parameters 
have been well constrained by the combined astronomical observations [4l|, • 

Subsequently, In [43j, Granda and Oliveros generalized the form of the IR cut-off on the basis of the Ricci scalar: 

L- 2 = aH 2 + /3H, (2) 
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where there are two independent model parameters a and j3, which can be determined by using the combined 
constraints of the thorough observational datasets. In this paper, we consider the holographic DE model with new 
IR cut-off in both fiat and non-flat case. The performance of a global fitting will be made by using the Markov Chain 
Monte Carlo (MCMC) method. In this way, we can work in the framework of multi-parameter freedoms, including 
the basic cosmological parameters (ft^h 2 , £l c h 2 , fifc) and the new-added model parameters (a, (3). 

The paper is organized as follows. In next section, we briefly review the holographic DE model with new IR cut-off. 
In section III, we perform the cosmic observation constraint on the holographic DE model. The last section is the 
conclusion. 



II. REVIEW OF HOLOGRAPHIC DARK ENERGY MODEL WITH NEW INFRARED CUT-OFF 

In this section, we give a brief review on the general formula in the holographic DE model with new IR cut-off. 
With a Friedmann-Robertson- Walker (FRW) metric 

rlr 2 

ds 2 = -dt 2 + a 2 {t)\-——r + r 2 {d9 2 + sin 2 9dtf)], (3) 
1 — fcr 

the Einstein field equation can be written as 

l_£(p 4 + 3P 4 ), (5) 



H 2 = 



a 

a = 6M 2 



where H is the Hubble function, and pi and Pi are the energy density and the pressure of a general piece of matter, 
and their subscripts i denote to, de and k, which respectively correspond to matter component, the holographic DE 
with new infrared cut-off and the curvature part of space. Here the matter component includes the cold dark matter 
and the baryon matter, i.e. 

Pm = Pcdm + Pb, P m = 0. (6) 

The parameter k = 1,0,-1 denote the closed, flat and open geometries, respectively. The effective energy density 
and the effective pressure of the curvature part are 

3M 2 k 

Pk = (7) 

Pk = ~Pk ~ ||. (8) 

As suggested by Granda and Oliveros in paper (43[, the energy density of the holographic DE with new IR cut-off 
is given as 

Pde = MI 2 pl {aH 2 +l3H), (9) 

where a and j3 are the dimensionless parameters in holographic DE model with new IR cut-off, which are regarded as 
independent of each other. In this paper, a dot denotes a derivative with respect to the cosmic time t. After changing 
the variable from the cosmic time t to x = In a, we can rewritten the Eq. (4) as 

ctMZi 2 ax 



With the help of the definitions as follows: 



ffo SM^H 2 ' H 2 ' K ' 



the Eq. (10) can be ulteriorly rewritten as 



1 dF 2 

E 2 = n k e~ 2x + n m0 e- 3x + aE 2 + -P—. (12) 

I dx 
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Solving this first order differential equation about E 2 , we can obtain 

7^2 o v-^x _i_ o v-3x i 2a — 3/3 3 a - p 2 2< °7 1) a: 

is = S2 fc e +i2 m0 e + — o — r ^ i2 «ioe +-5 —"fee + Joe 3 

3p — 2a + 2 p — a + 1 

= n k e- 2x + n mQ e- 3x + n de {x), (13) 

where /o is the integral constant and can be derived from the initial condition Eq = 1, which is / = 1 — /3 _^ +1 ^lfc — 
3ff-2«+2 ^ m0 ' anc ^ ^rfe(a;) is the dimensionless energy density of the holographic DE with new IR cut-off: 

^ , \ 2a - 3/3 ^ _o T a-/3 ^ _ 2r 1 2 2( a -i) 

"^ = 3/3-2a + 2 n ^ e + p^TT^ e + ( 1 -p^TT^- 3 / 3-2a + 2 ^ e ^ ' 

Then, combining the above definition of the dimensionless energy density of the holographic DE with its conservation 
equation, we can obtain the EOS of the holographic DE with new IR cut-off 

w de (z) = -l + — — . (15) 

In addition, we shall investigate the evolution of the deceleration parameter. Combining Eqs. (4) and (5) with the 
definition of the deceleration parameter, we can get 



2 2Y,iPi 
1 , 3 EiPi 



2 2n m0 {i + zf + n k {i + zf + n de { z y 



(16) 



where we have utilized the definitions of Pi — % and p~i = „ , „z . According to the energy conservation 

equation, we have 

5 ~ ld Pi ft^ 

Pi = - pi -3^- (17) 

In the review above, it is direct and natural to consider the parameters (3^0 and a ^ 1 when we solve the 
differential Eq. (12). Next, we discuss three special cases when the denominators in Eq. (13) equal zero as follows: 
Case 1: 3/3 - 2a + 2 = and (3 - a + 1 = 

In this case, we obtain (3 — and a = 1. Now the energy density of the holographic DE is p de — iM^H 2 . Compared 
to the Friedmann equation, it is found that the DE density is not consistent with the current observations of 70% 
exotic component in the universe. 

Case 2: 3/3 - 2a + 2 ^ and (3 - a + 1 = 

From the latter equation, we can get (3 = a — 1. Then the energy density of the holographic DE is p de = 
SM^i 9 ^-^— + aH 2 ). In this case, the solution of the first order differential equation is 



E 2 =f e~ 2x + JL-O m0e -3* - JL^ x n k e~ 2x 



=n k e~ 2x + n m0 e- ix + Q de (x), (18) 

where we have set the integral constant /o = Sl k and defined Q de (x) = 3 ^^Q m0 e~ 3x — -^r[xQ k e~ 2x as the dimensionless 
energy density of the holographic DE. Considering the current value of the dimensionless DE density, we can obtain 
the value of the only model parameter a = 3 ^ m0 _ff d d<i0 . So the case is a viable DE model. 
Case 3: 3/3 - 2a + 2 = and (3 - a + 1 ^ ff° 

From the former equation, we can obtain /3 = Then the energy density of the holographic DE is p de = 

3Mp Z (^-!-^|- + aH 2 ). In this case, the solution of the first order differential equation is 

E 2 =f e~ ix - -^—xn„ M e' 3x - -^—n k e- 2x 
a — 1 a — 1 

=n m0 e- 3x + n k e- 2x + n de ( x ), (19) 

where we have taken the integral constant f = £l m0 and defined Cl de (x) — — -^jxn m oe^ 3x — ^^l k e^ 2x as the 
dimensionless energy density of the holographic DE. Considering the current value of the dimensionless DE density, 
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TABLE I: The data fitting results of the cosmological parameters and the model parameters with la, 2a regions in the flat and 
non-flat holographic DE model with new infrared cut-off and ACDM model (flat case and non-flat case), where the combined 
observational data from SN 397, BAO and CMB and OHD are used. 



we can get the value of the only model parameter a = nj° n +n fc ■ So the case with the non-flat background geometry 
is a viable dark energy model. The result in the flat case is the same as Case 1. 

As far as the Case 2 and Case 3 discussed above are concerned, the model parameters a and j3 are reliant 
on each other. In this paper, we consider the generalized case with independent model parameters. In the next 
section, according to the combined observational data, we restrict the basic cosmological parameters and the model 
parameters, all of which are independent on each other. 

III. METHOD AND RESULTS 

In this section, we present the method and the data we have used. In our analysis, we perform a global fitting 
on determining the cosmological parameters using the Markov Chain Monte Carlo (MCMC) method. Since the 
computational requirements of MCMC procedures are insensitive to the dimensionality of the parameter space, we 
can expand the dimension of the parameter series, comparing with the traditional Maximum Likelihood (ML) method. 
The MCMC method is based on the publicly available CosmoMC package [44{, which has been modified to include 
the new parameters a and /? with having taken the weak priors as a e [0.5, 1.5] and (3 _ [0.1, 1.0]. Besides the two 
independent model parameters, the basic cosmological parameters are also varying with top-hat priors: the physical 
baryon density f^/i 2 _ [0.005,0.9], the dark matter energy density tt c h 2 £ [0.01,0.99], and in the non-flat case, the 
additional parameter flk € [—0.1,0.1]. In addition, we obtain three derived parameters f2„ e o, Q m o and the Hubble 
constant Hq from the basic cosmological parameters. 

In our calculations, we have taken the total likelihood L oc e _x I" 1 to be the product of the separate likelihoods of 
SN, BAO, CMB and OHD. Then the X 2 is 

X 2 = XSN + XBAO + XCMB + XOHD- (20) 

The expressions of x 2 s and datasets used in our paper are presented in Appendix lAl 

The best fit values of the cosmological parameters and the model parameters with la, 2a errors in holographic DE 
model with new IR cut-off and the ACDM model for the flat case and the non-flat case are listed in Table [I] We 
calculate the values of x 2 /dof, where dof is the compact notation of the number of degrees of freedom and equals 
the number of observational data points minus the number of free parameters. It is found that the values of x 2 /dof 
exhibit a significant difference between the flat case and the non-flat case in the holographic DE model with the new 
IR cut-off. In this two instances, it is seen that the non-flat holographic DE model with a smaller value of x 2 /dof is 
much supported by the current observations. Subsequently, comparing the value of x 2 /dof in the non-flat holographic 
DE model with those in the ACDM models, we find the differences are not obvious. From the minor differences among 
the values of x 2 /dof, we can conclude that the current combined datasets do not really favor the holographic DE 
model with the new IR cut-off over the concordance model. It is seen that the current observations support the flat 
concordance model with the smallest value of X 2 /dof the best. In Fig. [TJ [21 we show one dimensional probability 
distribution of each parameter and two dimensional plots for parameters between each other in the flat holographic 
DE model with new IR cut-off and the flat ACDM model. The corresponding plots in the non-flat holographic DE 
model with new IR cut-off and the non-flat ACDM model are presented in Fig. 13141 

Then, we investigate the evolutions of dark energy EOS and deceleration parameter in the holographic DE model. 
We consider the propagation of the errors for w(z) and q{z) by the Fisher matrix analysis. The errors are evaluated 
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FIG. 1: 1-D constraints on individual parameters (tlth 2 , a, f3, QdeO, SimO, -Ho) and 2-D contours on these parameters with la, 2a 
errors between each other using the combination of the observational data from SN 397, BAO, CMB and OHD in the flat 
holographic DE model with new IR cut-off. Dotted lines in the f-D plots show the mean likelihood of the samples and the 
solid lines are marginalized probabilities for the parameters in the flat holographic DE model with new IR cut-off [44] . 



by using the covariance matrix C« of the fitting parameters [4 a . |46[ , which is the inverse of the Fisher matrix and 
given by 

where 9 is a set of parameters, and lnL is the logarithmic likelihood function. The errors on a function / = f{9) in 
terms of the variables 9 are given by 0, |47j 



J 3 



where n is the number of parameters. Here, / will be dark energy EOS w{z;9i) or deceleration parameter q(z;9i). 
The parameters 9i respectively represent (il c h 2 , fife/i 2 , a, f3) for the flat case and (il c h 2 , fife/i 2 , flk, P) for the non-flat 
case. As shown in Fig. [S] (flat case) and Fig. [5] (non-flat case), we plot the evolutions of w(z) and q(z) with errors by 

w la {z) = w(z)\ e=s ± <j w , (23) 
qia(z) = q{z)\ e= g ±<J q , (24) 
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FIG. 2: 1-D constraints on individual parameters (Qbh 2 , £2ao, fimo, Ho) and 2-D contours on these parameters with la, la errors 
between each other using the combination of the observational data from SN 397, BAO, CMB and OHD in the flat ACDM 
model. Dotted lines in the 1-D plots show the mean likelihood of the samples and the solid lines are marginalized probabilities 
for the parameters in the flat ACDM model 



where are the best fit values of the constraint parameters. 

In Fig. [5] and Fig. [6j it is found that the combined observational data provide a fairly tight constraint on the 
holographic DE model with new IR cut-off. From the left panel in Fig. [5J it is seen that the EOS Wd e (z) with the 
best fit values can cross the boundary —1 and its current value is Wd e {z = 0) = —0.1414. From the right panel in Fig. 
[5j the value of the current deceleration parameter is given by qo = —0.7476 ± 0.0466. In Fig. [51 we get the current 
value of the dark energy EOS in the non-flat case is Wd e {z = 0) = —1.0653 < — 1, i.e. it is also phantom-like. The 
present value of the deceleration parameter is (70 = —0.6231 ± 0.0569. 



IV. CONCLUSION 



In summary, in this paper we have performed a global fitting on the parameters in the holographic DE model with 
new IR cut-off for the fiat case and the non-flat case, using a combined cosmic observations from type la supernovae, 
baryon acoustic oscillations, Cosmic Microwave Background and the observational Hubble data. The same constraints 
are performed on the flat and non-flat concordance models by using the same combined datasets. According to the 
Markov Chain Monte Carlo (MCMC) analysis, it is shown that the best fitting values of the model parameters (a, /3) 
in the flat holographic DE model with new IR cut-off tend to be smaller than those in the non-flat case. In the 
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FIG. 3: 1-D constraints on individual parameters {Q,bh 2 , Qk, «, fi, fideo, SimO, Ho) and 2-D contours on these parameters with 
la, 2a errors between each other using the combination of the observational data from SN 397, BAO, CMB and OHD in the 
non-flat holographic DE model with new IR cut-off. Dotted lines in the 1-D plots show the mean likelihood of the samples and 
the solid lines are marginalized probabilities for the parameters in the non-flat holographic DE model with new IR cut-off [44| ]. 



holographic DE models, the non-flat case with a smaller value of \ 2 /dof is much supported by the observations. 
In the non-flat cases, we have obtained the constraint values of the curvature terms Qk = O.OSOS^Q Q^io oi76 for 
the holographic DE model with new IR cut-off and Qk = — 0-0013^0 0076-^0 oios f° r ^ ae concordance model. These 
results indicate the two kinds of the non-flat background geometries in the two models. Then by using the best fit 
parameters, we plot the evolutions of the dark energy EOS and deceleration parameter with errors. From Fig. [5] 
and Fig. [51 it is found that the EOS of the holographic DE with new IR cut-off can cross the phantom divide —1, 
respectively with the current best values Wd e o = —1-1414 (fiat case) and Wdeo = —1.0653 (non-flat case). Comparing 
the flat and non-flat holographic DE models with the corresponding cases in the ACDM model, we can find that the 
current combined observations do not favor the holographic DE model with new IR cut-off over the ACDM model. 
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FIG. 4: 1-D constraints on individual parameters ({lf,h 2 , n*., JIao, QmO, Ho) and 2-D contours on these parameters with la, 2a 
errors between each other using the combination of the observational data from SN 397, BAO, CMB and OHD in the non-flat 
ACDM model. Dotted lines in the 1-D plots show the mean likelihood of the samples and the solid lines are marginalized 
probabilities for the parameters in the non-flat ACDM model [44| . 




FIG. 5: The evolutions of the EOS (the left panel) and the deceleration parameter (the right panel) in the flat holographic DE 
model with new IR cut-off with respect to the redshift z, where we have used the best fittings of the cosmological parameters 
and model parameters with la errors. 
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FIG. 6: The evolutions of the EOS (the left panel) and the deceleration parameter (the right panel) in the non-flat holographic 
DE model with new IR cut-off with respect to the redshift z, where we have used the best fittings of the cosmological parameters 
and model parameters with la errors. 



Appendix A: Cosmological Constraints Methods and Dataset 



1. Type la Supernovae constraints 

We use the SN la Constitution dataset, which includes 397 SN la (H- The 90 SN la from CfA3 sample with low 
redshifts are added to 307 SN la Union sample l2al. The CfA3 sample increases the number of the nearby SN la and 
reduces the statistical uncertainties. Following [48l l49j . one can obtain the corresponding constraint by fitting the 
distance modulus fj,(z) as 

Hth(z) = 51og 10 [D L (z)] + fiQ. (Al) 

In this expression Dl(z) is the Hubble-free luminosity distance Hodi,{z)/c, with Ho the Hubble constant, defined 
through the re-normalized quantity h as Hq = 100ft- km s - Mpc - 1 , and 

diW = W ston|vW I (A2) 

fio = 42.38 - 51og 10 /i, (A3) 

where smnn(-y/|fifc|x) respectively denotes sin(^/\Tlk\x) , y/\£lk\x, smh(y/\£lk\x) for < 0, = and fifc > 0. 
Additionally, the observed distance moduli fJ, bs(zi) of SN la at z, is 

Hobs{zi) = m obs (z l ) ~ M, (A4) 

where M is their absolute magnitudes. 

For the SN la dataset, the best fit values of the parameters p s can be determined by a likelihood analysis, based 
on the calculation of 

2/ HT'\ - \ ^ {Vobsjzi) - VthiPs, Zj)} 2 
X \Ps,M ) = 2^ 2 

SN Gi 

2 



^ i 5lo Sw[ D Up s ,z l )] - m bs(zj) + M'} 

SN 



where M' = fio + M is a nuisance parameter which includes the absolute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically }50i] as 



r+oa 

X 2 (Ps) = -2 In / exp 



— OO 



\x 2 {Ps,M') 



dM', 



to obtain 

B 2 , (C 



X 2 =A-^ + ln(-), (A6) 
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with 

, 2 



A _ {51ogi [-P£(p s ,^)] - m obs (zi)} 

^—^ n 2 



SN 1 

5 log 10 [D L {p s , Zi)] ~m obs {zi) 



SN 



SN "l 

Relation (|A5[) has a minimum at the nuisance parameter value M' — B/C, which contains information of the values 
of h and M. Therefore, one can extract the values of h and M provided one get the knowledge of one of them. Finally, 
it is noted that the expression 

X 2 SN (p s ,B/C) = A-(B*/C), 

which coincides to (|A6[) up to a constant, is often used in the likelihood analysis @ US [H|] , and thus in this case the 
results will not be affected by a flat M' distribution. 



2. Baryon Acoustic Oscillation constraints 



The Baryon Acoustic Oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, and measure the distance-redshift relation at z = 0.2. Additionally, Baryon Acoustic Oscillations in the 
clustering of the SDSS luminous red galaxies measure the distance-redshift relation at z — 0.35. The observed scale 
of the BAO calculated from these samples, as well as from the combined samples, are jointly analyzed using estimates 
of the correlated errors to constrain the form of the distance measure Dy(z) [271. l52ll53T| 



D v {z) 



(l + zfD\{z) 



cz 

W) 



1/3 



(A7) 



In this expression Da(z) is the proper (not comoving) angular diameter distance, which has the following relation 
with cIl(z) 



D a (z) 



d L (z) 
{l + zf 



(A8) 



The peak positions of the BAO depend on the ratio of Dv(z) to the sound horizon size at the drag epoch (where 
baryons were released from photons) z^, which can be obtained by using a fitting formula [54| : 



Zd 



1291(O m /i 2 



1 + 0.659{Q m h 2 )°- : 



-[i + h(n b h 2 ) 



with 



bi = 0.313(Q m /i 
b 2 = 0.238{n m h 



2 \ —0.419 I 



1 + 0.607{fl m h 



2n0.674i 



2\0.223 



(A9) 



(A10) 
(All) 



In this paper, we use the data of r a (zd)/ 'Dv(z) extracted from the Sloan Digitial Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) [Hj], which are listed in Table [TT1 where r s (z) is the comoving sound 
horizon size 



r s (z) =c / — — 



c 

75 Jo 



c s dt 
a 



a c,da 



a 2 H 



dz 



H{z) 



da 



a 2 H{a)^fl + (3n b /(4n 7 )a)' 
where c s is the sound speed of the photon— baryon fluid [55T - [57l | : 

4 



4 x Pb {z) 
3 p 7 (z) 



(A12) 



(A13) 
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z 


r s (z d )/D v (z) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE II: The observational r s (z d )/Dv(z) data [23] - 

and here fi 7 = 2.469 x Kr 5 /i~ 2 for T C mb = 2.725K. 

Using the data of BAO in Table llll and the inverse covariance matrix V^ 1 in [27j |: 

1 = ( 30124.1 -17226.9 \ 
^ -17226.9 86976.6 J ' 

thus, the Xbao(.Ps) 1S gi yen as 

X%Ao(Ps)=X t V- 1 X, (A15) 
where X is a column vector formed from the values of theory minus the corresponding observational data, with 

/ r s (z d ) _ _ 1905 33 \ 



X=[ °r±V J " , (A16) 



rs(zd) 
£>v(0.35) 



0.109715 



and X 1 denotes its transpose. 



3. Cosmic Microwave Background constraints 

The CMB shift parameter R is provided by [58| 



R(z*) = ^Jtt m H$(l ± z,)D A (z,)/c, (A17) 



which is related to the second distance ratio Da{z*)H(z*)/c by a factor y/X + z*. The redshift z* (the decoupling 
epoch of photons) is obtained using the fitting function [59j 

= 1048 [l + 0.00124(fi 6 /i 2 )- - 738 ] [1 + gi {n m h 2 ) 92 } , (A18) 

where the functions gi and g2 read 

gi = 0.0783(r! ;) /i 2 )- - 238 (l±39.5(r! 6 /i 2 ) - 763 ) _1 , (A19) 
52 = 0.560 (l± 21. l(fi b /i 2 ) 1 - 81 ) -1 . (A20) 

In additional, the acoustic scale is related to the first distance ratio, I?yi(z*)/r s (z*), and is defined as 

Z A = (1 + ^)Z^M. (A21) 

Using the data of R, I a, lOOilbh 2 and their covariance matrix of [-R(z»), Ia(z*), lOOf^/i 2 ] referring to [H,^, we can 
calculate the likelihood L as Xcmb = ~ 21nL: 

X 2 CMB = Ad^Cov-^d^dj)}^}', (A22) 

where Ad, ; = d. t — df ata is a row vector, and di = (R, I a, lOOf^/i 2 ). 
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4. Observational Hubble Data constraints 

The observational Hubble data are based on differential ages of the galaxies IrJOll. In [6l| , Jimenez et al. obtained an 
independent estimate for the Hubble parameter using the method developed in [601 ] . and used it to constrain the EOS 
of dark energy. The Hubble parameter depending on the differential ages as a function of redshift z can be written in 
the form of 



H(z) 



1 dz 



1 



dt 



(A23) 



So, once dz/dt is known, H(z) is obtained directly [64]. By using the differential ages of passively-evolving galaxies 
from the Gemini Deep Deep Survey (GDDS) [63[ and archival dat a [63-1691 . Simon et al. obtained H{z) in the range 
°f z ^ 1-8 [sH- The twelve observational Hubble data from |32l. |33| are list in Table Hill In addition, in (34|, 



z 





0.1 


0.17 0.27 


0.4 


0.48 


0.88 


0.9 


1.30 1.43 


1.53 


1.75 


H{z) (km s^Mpc -1 ) 


74.2 


69 


83 77 


95 


97 


90 


117 


168 177 


140 


202 


la uncertainty 


±3.6 


±12 


±8 ±14 


±17 


±60 


±40 


±23 


±17 ±18 


±14 


±40 



TABLE III: The observational H(z) data HHH. 

the authors took the BAO scale as a standard ruler in the radial direction, obtaining three more additional data: 
H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and H{z = 0.43) = 86.45 ± 3.27. 

The best fit values of the model parameters from observational Hubble data [62| are determined by minimizing 

XIM = t [g ^^~f^ )]2 , (A24) 

1=1 x ' 

where p s denotes the parameters contained in the model, H t h is the predicted value for the Hubble parameter, H f, s 
is the observed value, a(zi) is the standard deviation measurement uncertainty, and the summation is over the 15 
observational Hubble data points at redshifts 2^. 
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